0^ 
O 
O 
(N 



in 



o 



Probing neutron stars with gravitational waves 



Benjamin J. Owen 
The Pennsylvania State University 



LIGO-T0900053-V1 



Endorsed by: 
David H. Reitze (University of Florida) 
> : Stanley E. Whitcomb (LIGO-Caltech) 

O 

rn 
O 

a^ 
o 

> 



Introduction 



Within the next decade gravitational-wave (GW) observations by Advanced LIGO \i\ 
in the United States, Advanced Virgo [2] and GEO HF [32] in Europe, and possibly other 
ground-based instruments will provide unprecedented opportunities to look directly into the 
dense interiors of neutron stars which are opaque to all forms of electromagnetic (EM) ra- 
diation. The 10-10^ Hz frequency band available to these ground-based interferometers is 
inhabited by many neutron star mode frequencies, spin frequencies, and inverse dynami- 
cal timescales. GWs can provide information on bulk properties of neutron stars (masses, 
radii, locations. . . ) as well as microphysics of their substance (crystalline structure, viscos- 
ity, composition. . . ), some of which is difficult or impossible to obtain by EM observations 
alone. The former will tell us about the astrophysics of neutron stars, and the latter will 
illuminate fundamental issues in nuclear and particle physics and the physics of extremely 
condensed matter. Although GW searches can be done "blind," they become richer and 
more informative with input from EM observations; and thus the combination of the two 
is crucial for learning the most we can about neutron stars. Healthy GW and EM observa- 
tional programs must be accompanied by vigorous theoretical research on the interface of 
astrophysics, gravitational physics, nuclear and particle physics in order to extract the most 
from the observations. 

Advanced LIGO has been funded and it and its international partners will not be ranked 
as part of the decadal review. Our purpose here is rather to describe the neutron-star 
aspect of GW science and point out observational and theory issues in EM astronomy and 
astrophysics which connect to it, and thus are important to maximizing the scientific output 
of the advanced detectors (LIGO, Virgo, and GEO HF) and the future underground GW 
detectors now in planning. 

Opportunities to learn about neutron stars through GW and and GW/EM observations 
naturally divide into four categories: supernovae, binary mergers, starquakes, and continuous 
waves. The first is addressed in another white paper |13j; here we address the other three. 

Binary Mergers 

Advanced LIGO is likely to observe mergers of double neutron star (NS/NS) binaries at 
a rate of a few to a few hundred per year; and black-hole/neutron-star (BH/NS) binaries 
perhaps in a comparable range of rates [18]. If the observed rates are at either extreme 
of the range, they would refiect on extremes of compact binary populations and aspects of 
their formation such as the properties of a common envelope phase and the distribution of 
neutron-star birth kicks, both of which can prematurely merge or disrupt binaries before the 
second neutron star or the black hole is formed. 

Individual merger signals can reveal aspects of neutron star structure: For BH/NS merg- 
ers the tidal disruption of the neutron star can take place at frequencies where ground- 
based GW detectors are most sensitive. The latest, fully relativistic, numerical simulations 
of merger [TJ] confirm the earlier Newtonian intuition that the "break frequency" of the 
merger signal, the frequency at which much of the neutron star is swallowed by the black 
hole and the signal becomes very faint, is strongly correlated with the neutron star radius 
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and therefore can serve as a good measurement of it. Beyond Advanced LIGO, or with a 
very lucky event during the time of Advanced LIGO, a neutron star's mass can be measured 
to as well as 20% from the precessional modulation of the signal if the companion is a rapidly 
rotating black hole [27j; and a high neutron star mass would significantly constrain the star's 
composition. Recent analytical [16] and numerical estimates |2lj of the effects of tidal defor- 
mations on NS /NS orbits (and thus the phase evolution of merger signals) indicate that the 
strongest signals observed by Advanced LIGO might yield measurements of stellar radii to a 
precision as good as 1 km, even better than the break frequency. Such a radius measurement 
constrains to order 10% the pressure around twice nuclear density, a number unattainable 
in terrestrial laboratories which is correlated with properties such as the incompressibility 
and isospin asymmetry energy of nuclear matter. 

GW and EM observations of mergers are important to each other as described in the 
transient white paper [13]: The availability of EM triggers localized in time and sky position 
increases the sensitivity of GW searches, and comparison of GW and EM observations can 
constrain properties of the merger and gamma-ray burst (if any). EM observations are im- 
portant to GW mergers in another respect: While there are several observed NS/NS systems 
on which to empirically base merger rates, and eventually make population statements by 
comparing to observed GW merger rates, the BH/NS event rates are estimated purely based 
on theoretical population synthesis codes which contain many poorly known factors. If a 
pulsar survey, performed for instance by the Square Kilometer Array [3], reveals any BH/NS 
binaries it would be of great relevance to GW merger observations. 

In order to get the most science out of the observations it is also crucial to support 
a program of numerical simulations of increasing complexity and realism. Details of the 
neutron star matter such as equation of state, viscosity and other transport coefficients, 
photon and neutrino luminosities are modeled at present very simply; yet they may have 
important signatures in the GW emission and certainly are crucial to the EM emission. It 
is also important to continue work on analytical and semi-analytical perturbation theory, as 
witnessed by the recent analytical prediction of an important tidal effect [TB] which preceded 
numerical investigations ^24j due to the computational cost of lengthy simulations. Advances 
in population simulation [18] are also important to extract astrophysical information from 
GW event rates. 

Star quakes 

LIGO has already placed upper limits on GW bursts associated with flares of highly 
magnetized neutron stars or magnetars [5| i9j, in some cases comparable to the 10^^ erg 
associated with the largest giant flare; and Advanced LIGO will be sensitive to GW energies 
more than two orders of magnitude lower. The dominant model of the flares [25] associates 
them with starquakes due to sudden rearrangements of the neutron star's intense magnetic 
fleld, quakes which would undoubtedly excite the star's quasinormal modes and thus emit 
GWs. Indeed the quasiperiodic oscillations seen in the x-ray tails of giant flares may be 
modes of the crust, influenced by the magnetic fleld [17j. The GW upper limits are already 
interesting in the sense that they are just entering the range predicted by the most extreme 
theoretical models. However, it is not clear what range of GW energies would be associated 
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with other models of the flare, because there has been relatively little theoretical work on 
the subject. Thus there is a need for modeling to extract more meaning from current GW 
observations, let alone those in the next decade. 

GW searches so far have focused on frequencies tied to the quasiperiodic oscillations seen 
in x-rays after the flares, on ringdowns of the fundamental or /-modes which are probably 
the most efficient GW emitters, and opportunistic searches for random bursts around 100 Hz 
where the detectors' noise is low. The frequency of a detected /-mode would strongly con- 
strain the structure of the star — neutron stars would be easily distinguishable from quark 
stars, and the mean density can be measured to a percent |20]. The total energy radiated 
would reveal the mechanical efficiency of the flare mechanism and have something to say 
about how much of the action happens inside the neutron star as opposed to the magneto- 
sphere surrounding it. The energy and rate at which energy is transferred to the /-mode 
from the initial excitation could also reveal something about how the Alfven continuum of 
magnetic modes is coupled to the mechanical modes, and thus how the magnetic fleld is 
distributed within the star; and any timing offsets between GW and EM events would also 
constrain the possible couplings. Even more exciting would be a detection of a random 
burst not associated with a known mode or x-ray quasiperiodic oscillation: Such a signal 
could be evidence of turmoil of the magnetic fleld interior to the star rather than in the 
magnetosphere. 

The story will be similar for searches for GW bursts associated with pulsar glitches. Al- 
though the glitch mechanism still remains a puzzle, it can be reasonably expected that the 
transfer of angular momentum to the crust from the core, which is believed to be differ- 
entially rotating with respect to the crust, would excite quasinormal modes and thus GW 
emission. Extrapolating from [H]: With advanced ground-based interferometers tuned to 
kHz frequencies, modes excited by glitches of the nearby Vela pulsar might be observable; 
and future underground detectors would be sensitive to large glitches from pulsars within a 
few kpc. 

For both starquake scenarios as for mergers, EM astronomy plays a crucial role increasing 
the GW search sensitivity by localizing the times and sky positions of the events. Burst 
searches using non-matched flltering techniques also improve when provided with frequency 
windows, such as determined by observations of quasiperiodic oscillations or predictions of 
quasinormal mode frequencies. Therefore in order to help GW observations it is important to 
sustain x-ray, gamma-ray, and radio timing and rapid response capabilities. And to extract 
physical information on the interiors of neutron stars it is important to have more modeling 
of GW and EM emission mechanisms, and the phenomenology of how they relate to each 
other. 

Continuous Waves 

Long-lived continuous- wave GW signals may be the hardest to detect, and certainly they 
are the most computationally costly to search for; but with 10^-10^° cycles per year they 
could ultimately yield the most precise information on neutron stars. Here the connec- 
tions between GW and EM astronomy are particularly tight: Because all continuous- wave 
searches except monitoring of known (EM) pulsars are severely computationally limited. 
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electromagnetic observations can greatly increase the sensitivities of GW searches. The dif- 
ference between the best upper limits for all-sky and known-pulsar searches on the same 
set of LIGO data is nearly an order of magnitude in strain amplitude pi, [7] or a factor 100 
difference in luminosity. The all-sky survey also used much more computing power to achieve 
its lesser sensitivity. 

All-sky GW surveys can reveal entirely new neutron star populations. Of the 10*^- 
10^ neutron stars formed in our galaxy over 10^° years, only 2000 have been identified 
electromagnetically — the vast majority as pulsars — and only 20,000 could be identified as 
pulsars even with the Square Kilometer Array This is because most neutron stars are 
too old to pulse and pulsar EM emission is relatively narrowly beamed, meaning that only a 
fraction of currently EM-active pulsars can be observed. But GW emission is beamed very 
little, and the emission mechanisms (which ultimately rely on some nonaxisymmetry of the 
mass or momentum distribution in the star) may not be limited in activity to the youth 
of the star or a possible "recycled by accretion" phase as they are for pulsars. The main 
selection effect is that continuous GW emission is detectable only for spin frequencies in the 
ground-based detector band, about 1/5 of known pulsars for Advanced LIGO and Advanced 
Virgo. Although population studies of continuous wave sources are still in their infancy [19j, 
LIGO searches are already coming within striking distance of the most extreme scenarios 
and Advanced LIGO will have something to say about more likely scenarios. 

Coherent integration of a year of data can allow LIGO to localize a continuous wave 
source's sky location to the sub-arcsecond level [14], facilitating followup observations to 
find counterparts in x-rays, optical, and other electromagnetic wavebands. Gravitational- 
wave timing of newly detected spinning neutron stars will also aid followup searches for 
electromagnetic pulsations by providing spin frequencies and frequency derivatives. An ex- 
tended observation with Advanced LIGO can also, through sphericality of the wavefronts, 
obtain a distance measurement (independent of EM observation) to about 10% [25] . 

A distance measurement (whether GW-based or EM-based from a known pulsar) allows 
a determination of the neutron star's quadrupolar deformation or ellipticity, a quantity of 
great interest since it may shed light on whether the neutron star is indeed made of neutrons 
or contains quarks or other exotica. Very large ellipticities are sustainable in quark models 
but not in normal nuclear matter, and thus detection of a large enough ellipticity would 
confirm not only the existence but also the crystalline nature of quark matter [22] . 

Rapidly accreting neutron stars (whether pulsing or not) in low-mass x-ray binaries are 
also very interesting targets for continuous GW searches because the accretion is known 
to be asymmetric due to the occasional observation of x-ray pulsations in some systems. 
This asymmetry can lead to a GW-emitting mass asymmetry through several mechanisms 
including the temperature-dependent electron capture onto nuclei in the crust [12], mag- 
netic funneling of accreted material [21], and sustained instability of rotational "r-mode" 
oscillations in the fluid below the crust [TT] (which are likely to have a briefer episode of 
continuous GW emission in young neutron stars [23j). The observed lack of extremely high 
spin frequencies in the most rapidly accreting neutron stars also strongly suggests that their 
accretion torques are balanced by GW emission torque [12]. In the simplest version of this 
torque balance scenario, Sco X-1 would be detectable by advanced interferometers, while 
many systems would be detectable with future underground instruments; and if the GW 
emission (like the x-ray emission) also goes through active and quiescent phases, many sys- 
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terns would be detectable by advanced interferometers if they go through a GW-active phase 
in the next decade ^Uj. Determination of the neutron star spin period in Sco X-1 or other 
systems where it is completely unknown increases the sensitivity of a GW search by a factor 
of a few in strain amplitude, or an order of magnitude in luminosity. Detection of a signal 
from an accreting neutron star would confirm the hypothesis that the observed spins of older 
neutron stars are due to GW/accretion torque balance, and would constrain the mechanics 
of accretion as well. 

Particularly interesting is the case when GW and EM observations are made simultane- 
ously, whether in pulsars or in the many accreting neutron stars where the spin frequency is 
only roughly known (for instance through x-ray burst oscillations): The ratio of gravitational- 
wave frequency to spin frequency immediately identifies the emission mechanism, whether 
a static deformation (ellipticity) of the rotating star or an oscillating r-mode. In the latter 
case, the frequency ratio has some dependence on the neutron star structure [5] and thus, 
once observed, can be used to measure the equation of state of dense matter. Due to a com- 
plicated interplay of factors involving the microphysics of the matter deep inside the star 
and possibly its effect on fluid and superfluid dynamics, detection of a long-lasting r-mode 
signal from an accreting neutron star (at any frequency) is a strong indicator for the pres- 
ence of some type of strange matter (quarks, mesons, or hyperons) in the core (e.g. [23 [TU]). 
If precise EM timing is available on long timescales, even more is possible: A comparison 
between GW and EM timing can constrain the coupling between the solid GW-emitting 
component of the star and whatever component of the star and magnetosphere is emitting 
EM [8]; and this comparison across a glitch may shed light on the glitch mechanism. Future 
underground interferometers might also be able to observe (for the first few days after a 
glitch) a continuous GW signal associated with the fluid interior's response to the abrupt 
relative motion of the crust; and this signal would carry information on the superfluidity 
and transport coefficients of the matter in the interior [28| . 




What is needed to realize this potential? Since searches for GW from known pulsars rely 
on coherent phase models lasting of order a year, it is crucial to have regularly updated EM 
pulsar timing data synchronous with LIGO and Virgo data runs. With most pulsars pulsing 
in radio, it is important to maintain healthy radio observatories and monitor those pulsars 
with high spindowns which leave room for strong GW emission — such as PSRs J1952+3252 
and J1913-I-1011, which have not been observed for some time. Also one of the most rapidly 
spinning down and frequently glitching pulsars, J0537-6910, can only be timed at the mo- 
ment with the Rossi X-ray Timing Explorer satellite. The extreme glitchiness of this pulsar 
indicates a substantial and quickly changing solid component, which also indicates likely 
strong gravitational-wave emission. However it also means that without regular timing a 
LIGO search for J0537 would quickly become intractable, as the number of possible timings 
to be used in the analysis would grow out of control. Since RXTE is scheduled to go down 
soon, it is important to get a comparable mission back up quickly. It is also helpful to GW 
observations to discover as many new pulsars as possible (in radio, x-rays, or any band), 
and to discover pulsations in stars where they have not yet been observed. Further hunts for 
young neutron stars and massive star forming regions will also help GW searches by produc- 
ing more targets for point or small-area searches, which can be more sensitive than all-sky 
searches [31J- Finally, it is crucial to encourage work on the theory of emission mechanisms 
and all the complicated physics of neutron stars, as this is needed to turn the combined EM 
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and GW data streams into statements about the properties of matter in states unreachable 
in terrestrial laboratories. 



References 

[1] http://www.ligo.org/. 

[2] http://www.virgo.infin.it/. 

[3] http://www.skatelescope.org/. 

[4] B. Abbott et al. All-sky LIGO Search for Periodic Gravitational Waves in the Early S5 
Data. arXiv:0810.0283. 

[5] B. Abbott ct al. Search for gravitational wave radiation associated with the pulsating 
tail of the SGR 1806-20 hyperflare of 27 December 2004 using LIGO. Phys. Rev., 
D76:062003, 2007. 

[6] B. Abbott et al. Upper limits on gravitational wave emission from 78 radio pulsars. 
Phys. Rev., D76:042001, 2007. 

[7] B. Abbott et al. All-sky search for periodic gravitational waves in LIGO S4 data. Phys. 
Rev., D77:022001, 2008. 

[8] B. Abbott et al. Beating the spin-down limit on gravitational wave emission from the 
Crab pulsar. Astrophys. J., 683:L45-L50, 2008. 

[9] B. Abbott et al. Search for Gravitational Wave Bursts from Soft Gamma Repeaters. 
Phys. Rev. Lett, 101:211102, 2008. 

[10] N. Andersson, D. I. Jones, and K. D. Kokkotas. Strange stars as persistent sources of 
gravitational waves. Mon. Not. Roy. Astron. Soc, 337:1224, 2002. 

[11] N. Andersson, K. D. Kokkotas, and N. Stergioulas. On the relevance of the r-mode 
instability for accreting neutron stars and white dwarfs. Astrophys. J., 516:307, 1999. 

[12] L. Bildstcn. Gravitational radiation and rotation of accreting neutron stars. Astrophys. 
J., 501:L89, 1998. 

[13] J. S. Bloom et al. Astro2010 Decadal Survey Whitepaper: Coordinated Science in the 
Gravitational and Electromagnetic Skies. arXiv:0902.1527. 

[14] P. R. Brady, T. Creighton, C. Cutler, and B. F. Schutz. Searching for periodic sources 
with LIGO. Phys. Rev., D57:2101-2116, 1998. 

[15] Z. B. Etienne et al. Fully General Relativistic Simulations of Black Hole- Neutron Star 
Mergers. Phys. Rev., D77:084002, 2008. 



6 



[16] E. E. Flanagan and T. Hinderer. Constraining neutron star tidal Love numbers with 
gravitational wave detectors. Phys. Rev., D77:021502, 2008. 

[17] G. Israel et al. Discovery of Rapid X-ray Oscillations in the Tail of the SGR 1806-20 
Hyperfiare. Astrophys. J., 628:L53-L56, 2005. 

[18] V. Kalogera, K. Belczynski, C. Kim, R. O'Shaughnessy, and B. Willems. Formation of 
Double Compact Objects. Phys. Rept, 442:75-108, 2007. 

[19] B. Knispel and B. Allen. Blandford's Argument: The Strongest Continuous Gravita- 
tional Wave Signal. Phys. Rev., D78:044031, 2008. 

[20] K. D. Kokkotas, T. A. Apostolatos, and N. Andersson. The inverse problem for pulsating 
neutron stars: A 'fingerprint analysis' for the supranuclear equation of state. Mon. Not. 
Roy. Astron. Soc, 320:307-315, 2001. 

[21] A. Melatos and D. J. B. Payne. Gravitational radiation from an accreting millisecond 
pulsar with a magnetically confined mountain. Astrophys. J., 623:1044-1050, 2005. 

[22] B. J. Owen. Maximum elastic deformations of compact stars with exotic equations of 
state. Phys. Rev. Lett, 95:211101, 2005. 

[23] B. J. Owen et al. Gravitational waves from hot young rapidly rotating neutron stars. 
Phys. Rev., D58:084020, 1998. 

[24] J. S. Read et al. Measuring the neutron star equation of state with gravitational wave 
observations. arXiv:0901.3258. 

[25] N. Seto. Gravitational Wave Astrometry for Rapidly Rotating Neutron Stars and Esti- 
mation of Their Distances. Phys. Rev., D71:123002, 2005. 

[26] C. Thompson and R. C. Duncan. The soft gamma repeaters as very strongly magnetized 
neutron stars - I. Radiative mechanism for outbursts. MNRAS, 275:255-300, July 1995. 

[27] M. V. van der Sluys ct al. Gravitational- Wave Astronomy with Inspiral Signals of 
Spinning Compact- Object Binaries. Astrophys. J. Lett., 688:L61-L64, 2008. 

[28] C. A. van Eysden and A. Melatos. Gravitational radiation from pulsar glitches. Class. 
Quant. Grav., 25:225020, 2008. 

[29] R. V. Wagoner. Conditions for steady gravitational radiation from accreting neutron 
stars. Astrophys. J., 578:L63-L66, 2002. 

[30] A. Watts, B. Krishnan, L. Bildsten, and B. F. Schutz. Detecting gravitational wave 
emission from the known accreting neutron stars. Mon. Not. Roy. Astron. Soc., 389:839- 
868, 2008. 

[31] K. Wcttc ct al. Searching for gravitational waves from Cassiopeia A with LIGO. Class. 
Quant. Grav., 25:235011, 2008. 

[32] B. Willke et al. The GEO-HF project. Class. Quant. Grav., 23:S207-S214, 2006. 



7 



